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ABSTRACT 

We investigate the scattering and absorption of hght by random baUistic aggre- 
gates of spherical monomers. We present general measures for the size, shape, and 
porosity of an irregular particle. Three different classes of ballistic aggregates are con- 
sidered, with different degrees of porosity. Scattering and absorption cross sections are 
calculated, using the discrete dipole approximation (DDA), for grains of three compo- 
sitions (50% silicate and 50% graphite; 50% silicate and 50% amorphous carbon; and 
100% silicate, where percentages are by volume), for wavelengths from 0.1 /im to 4/im. 
For fixed amount of solid material, increased porosity increases the extinction at short 
wavelengths, but decreases the extinction at wavelengths long compared to the overall 
aggregate size. Scattering and absorption cross sections are insensitive to monomer size 
as long as the constituent monomers are small compared with the incident wavelength. 
We compare our accurate DDA results with two other approximations: the analytical 
multi-layer sphere (MLS) model and effective medium theory (EMT). For high porosity 
and/or absorptive materials, the MLS model does not provide a good approximation for 
scattering and absorption by ballistic aggregates. The EMT method provides a much 
better approximation than the MLS model for these aggregates, with a typical differ- 
ence < 20% in extinction and scattering cross sections compared with DDA results, for 
all types, compositions and wavelengths probed in this study. 

Subject headings: dust, extinction - scattering - circumstellar matter - interplanetary 
medium 

1. Introduction 

The appearance of star- forming galaxies is strongly affected by interstellar dust. Starlight is 
absorbed and scattered by dust particles, and the absorbed energy is predominantly reradiated at 
infrared wavelengths. The optical properties of the dust must be characterized in order to determine 
the spectra and spatial distribution of the stars (i.e., to "correct" for reddening and extinction), 
and to interpret the observed infrared emission. Characterizing the dust is also important for more 
fundamental reasons: dust grains play an important role in the thermodynamics and chemistry of 
the interstellar medium (e.g., heating from photoelectric emission, and grain catalysis of H2), in 
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interstellar gas dynamics (e.g., radiation pressure on dust grains), and in the formation of stars, 
planets, and planetesimals. 

Many interplanetary dust particles (IDPs) are irregular, fluffy aggregates, sometimes described 
as havin g a "fractal" appearance. Direct evi dence comes from such IDPs collected by high-flying 
aircraft (jBrownled Il985l : IWarren et al.lll994l ). Laboratory and microgravity experiments of dust 
particle interactions, which may mimic the conditions in the early solar system, also suggest that 
particles form fractal assemb lies through ballistic aggreg ation dWurm Blumlll998l :l I Blum &: WurmI 
200d : iKrause &: Blumll2004l ). Similar aggregation processes may well occur in interstellar environ - 
ments, where small dust grains coagulate in dense molecular clouds (jDorschner &: Henningj|l995l ). 
A number of authors have proposed that interst ellar grains consist priin arily of such aggregate 
structures, with a random, porous geometry (e.g.. iMathis &: Whiffenlll989l ). 



Aggregates are irregular, porous, particle assemblies, possibly incorporating multiple materials. 
Early attempts to estimate the optical properties of such structures used a combination of Mie 
theory and effective medium theory (EMT) , where an aggregate is approximated by a homo geneous 
sphere with an "effective" r efractive index i ntended to allow for the e f fects of porosity (e.g.. iMathis 
19961 : iLi &: Greenberg||l998l ). Alternatively, IVoshchinnikov &: Mathid (|l999l ) proposed a multi-layer 
sphere (MLS) model to account for porous, composite dust grains. The advantage of this MLS 
model is that, lust as for hoin ogeneous spheres, the solution is exact and can be obtained via a 
fast algorithm (IWu et al.l 119971 ). The limitations are: the grains are assumed to be spherical, and 



all compositions (including vacuum) are assume d to be well mixed inside the sph ere. Comparisons 
of the MLS with EMT- Mie theory are done by IVoshchinnikov et al.l (120051 . l2006l ) , suggesting that 
the MLS model is accurate for spherical, well-mixed dust grains, for certain compositions and 
porosities. So far there has been no study of whether or not the MLS prescription is also suitable 
for irregular aggregates. 

Solving Maxwell's equations becomes increasingly challenging when the geometry departs from 
spheres. Several methods have been developed to compute light scattering by non-spherical parti- 
cles. The "ex tended boundary condition meth od" (EBCM) introduc ed by IWatermanI (jl97ll ). and 
developed by iMishchenko &: Travid (119941 ) and IWielaard et al.l (119971 ) can be applied to solid tar- 
gets with relatively smooth surfaces, and exact s eries expansions have been developed for spheroids 



(Asano &: Yamamotdll975l: lAsano &: Sato 



198C 



Voshchinnikov &: Farafonov 



1993 



Farafonov et al 



19961 ). but these techniques cannot be applied to the complex geometries of interest here. The 
special case of clusters o f spheres c an also be treated by superposition of vector spherical har- 



monics ( Mackowski 



1991 



Xul 119971 ) , sometimes referred to as the generalized multisphere Mie 



(GMM) solution. If the complete T-matrix can be found for the individual particles, the T-matrix 
for a cluster can, in principle, be found using the superposition T-matrix method (TMM, e.g., 
Mackowski &: Mishchenkd Il996l ). However, these techniques, although formally exact, can prove 
very computationally-demanding when applied to spheres that are numerous and in close proximity 
(e.g., in contact), particularly when the refractive index differs appreciably from vacuum. 
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The discrete dipole approximation (DDA) method (e.g., Purcell &: Pennvpackei|[l973 : Draine &: Flatau 



1994 ) is the most flexible among all the methods, although substantial computational resources 
may be needed in order to achieve the desired accuracy. Using either GMM, TMM or DDA, 
the optical properties of difi'erent kinds of aggrega,tes have been investigated extensively d uring 
the past decade (e .g., IWestlll99ll : iLumme k. Raholalll994l : IPetrova et al.ll200d : iKimura et al.ll2006i : 
Bertini et al.l 120071 ): most of these investigations are applied to cometary dust. 



There is increasing interest in the possibility that interstellar and circumstellar grains may be 
porous, random aggregates. To confront this hypothesis with observational data, it is necessary to 
calculate the optical properties of such aggregates. The calculations must be carried out for a wide 
range of wavelengths - from infrared to vacuum ultraviolet - and for aggregates with potentially 
realistic geometry and composition. Massive computations are required to attain these goals. 

In this paper, we use the DDA to calculate scattering and absorption for a large sample of 
aggregates with different compositions, sizes, and porosities, over a wide range of wavelength. We 
explore the dependence of various optical properties on wavelength, grain size, aggregate geometry, 
as well as compositions, and we investigate the applicability of multi-layer sphere calculations and 
the EMT-Mie model to these random aggregates. 

Characterization of the size, gross shape, and porosity of an irregular particle is discussed in 
f|2j ^describes procedures for generating random aggregates ("BAMl" and "BAM2") that are 
less porous than the standard ballistic aggregates, and ^H] discusses the applicability of the DDA 
method to these aggregates. 

Extinction and absorption behaviors are presented in ^ with comparison to the MLS model in 
§5.31 s-iicl comparison to the EMT-Mie model in §5.41 Our results are summarized and discussed in 
^ Ap plications to circ umstellar debris disks and cometary dust will be presented in a companion 
paper (IShen et al.l 120071 . hereafter Paper II). 



2. Geometric Properties of Irregular Particles 



In general, we define the "effective radius" (or "volume equivalent radius") of a general 3- 
dimensional structure to be the radius of a sphere with volume equal to the volume V of solid 
material in the structure: 

An 



1/3 



(1) 



There are several possible ways to quantify the "porosity" or openness of the structure, de- 
pending on how to measure its apparent size. Different definitions of the apparent size f or a fiuffy 
structure have been used, i.e., b ased on the radi us of gyration (e.g., iKozasa et al.lll992l ) or using 
the geometric cross section (e.g.. IOssenkopflll993l ). 



Here we propose a simple but general way to quantify the size, shape and porosity of a finite 
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structure with arbitrary geometry. Let the "density" pQ = 1 at points within sohd material, 
and pq = otherwise!^ Let I be the moment of inertia tensor of the structure, with eigenvalues 
Ii ^ I2 ^ H- It is convenient to define dimensionless quantities 



a,; 



0.4Ma2g 



(2) 



where M = pqV . For a target that is a solid sphere, = 1. 



We characterize the size and shape of an irregular structure by considering an "equivalent 
ellipsoid" of uniform density /po and semi-major axes a < b < c. The quantities (/, a, b, c) are 
uniquely determined by requiring that the mass M and principal components of the moment of 
inertia /j of the target and its equivalent ellipsoid be identical, such that: 



Sir . 
— PoOefr«i 

Svr . 
— Poaefr«2 

Stt 5 
-poa^gai 



with solutions 



a 
b 
c 
f 



15' 



(02 + as — ai 



fpo^abc{b^ + c^) 
15 

fpo—abc{a^ + c^) 
15 

fpo—abc{a^ + b^) 
15 



|l/2 
1/2, 



«eff 



(as + ai - a2) ' aeff 
(ai + a2 - as)^''^aeff 

[(a2 + as — ai)(as + ai — a2)(ai + a2 — as)]^"*^^^ 



(3) 
(4) 
(5) 
(6) 



(7) 
(8) 
(9) 
(10) 



/ from eq. (|lUp is then the volume filling factor for the equivalent ellipsoid. We define a characteristic 



size 



Rabc = {abcf/^ = [(a2 + as - ai)(as + ai - a2)(ai + a2 - as)]'^" ffloff (H) 

which is simply the radius of a sphere with a volume equal to the equivalent ellipsoid defined above; 
the filling factor is just / = {aes/Rabc)^ ■ We define the "porosity" to be 

V = l-f = l- ia,s/Rabcf . 



11/6, 



(12) 



Kozasa et all (j 19921 ) characterized the overall size of the structure by a radius proportional to the 
radius of gyration: 

"ai + a2 + as 



Rkbm — 



.1/2 



Oeff 



(13) 



We are discussing purely geometric properties of the target, which are independent of the actual material density. 
The function po is here introduced to allow use of familiar concepts like total mass and moment-of-inertia tensor, but 
should be understood to indicate "occupation" if po — 1, and vacuum if po = 0. 
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For a uniform density target, i?KBM = (5/3)^/^i?gyr, where Rgy^ is the radius of gyration. iKozasa et al 
(|l992l ) proposed that the porosity be based on the ratio between i?KBM and aes'- 

^kBM = 1 - (Oeff/^KBM)^ (14) 

Finally, one can also define a radius 

Rproj = [(Cproj)/^]'/' , (15) 

where (Cproj) is the orientation-averaged projected area of the target. For irregular targets, i?proj 
is more difficult to compute than either Rate or -Rkbm- These three radii - Rate, Rkbm, and -Rproj 
- will be compared below for the ballistic aggregates in this study. 

The gross shape of the target can be characterized by the ratios c/b and b/a, where a, b, c are 
given by eq. ([71l9|). We will refer to agglomerates as oblate if c/b < b/a, and prolate if c/b > b/a. 

It is convenient to define dimensionless efficiency factors Q for absorption, scattering and 
extinction: 

Q(A)^^ , (16) 

where C(A) is the cross-section for absorption, scattering, or extinction, and a^s is the solid- volume- 
equivalent radius, defined by eq. ([I]). 



3. Geometric Properties of Ballistic Aggregates 



We use 3 simple algorithms for generating irregular porous structures with varying degrees 
of porosity. The first is standard "ballistic agglomeration" (BA) , sometimes referred to as "bal- 
listic particle-cluster agg l omer a tion" (BPCA), previously discussed by many authors ( e.g., [West 
199ll : iKozasa et al]|l992l . Il993l : lOssenkop j ^993 : Iximura et aPbood : iBertini et alJboOTl l. We also 
introduce two new prescriptions for agglomeration - BAMl and BAM2 - that produce random 
structures that are less "fluffy" than those produced by BA (see ^S.ip . 



Random aggregates can be characterized by the Haussdorf dimension D (sometimes referred 
to as the fractal dimensionality) , with mass M oc R^ as i? — > cxd where R is some characteristic size 
of the agglomerate. The filling factor / cx M/R^ oc M^~^/^ . All three agglomeration algorithms 
employed here (BA, BAMl, BAM2) are thought to have Haussdorf dimension Z) « 3, with the 
filling factor / constant as M ^ cxo. 

Other processes have also been proposed for generating random agglomerates, including ballis- 
tic cluster-cluster aggregation" (BCCA) and diffusion-limited aggregation (DLA). BCCA clusters 
obtained by random , aggregation of clusters with equal-mass clusters are very "fluffy" , with D ~ 2.25 
(jKozasa et al.lll993l ) and filling factor / oc M"'^-^^. Clusters produced by diffusion- limited aggre- 
gation have a dendritic appearance, with D ~ 2.5 (jWitten &: Gated 119861 ) . and / oc M~^''^. We 
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do not consider either BCCA clusters or DLA clusters here, because their fragile geometries seem 
unlikely to be representative of circumstellar or interstellar grains, which are subject to occasional 
grain-grain collisions. 

We stress that we do not claim that any of our algorithms (BA, BAMl, or BAM2) provides 
a realistic representation of the actual processes responsible for the growth of circumstellar or 
interstellar grains - they are merely convenient procedures for generating irregular targets that 
may bear some geometric resemblance to real circumstellar or interstellar grains. 



3.1. Target Generation 



We construct clusters ("targets") by random ballistic agglomeration. Each cluster is composed 
of a certain number of spheres (or "monomers") of a single radius ao- The agglomeration process 
begins with one monomer (j = 1). A cluster is built up by sequential arrival of additional monomers 
(j = 2, ...,A^) on random rectilinear trajectories. We consider three different classes of clusters, 
distinguished by aggregation rules. 

BA ("ballistic agglomeration") clusters are produced by requiring arriving monomers to "stick" 
at th e poirit where they first contact the preexisting aggregate. This is a well-established procedur e 



(e.g.. I West 



1991 



Kozasa et al.lll992l . Il993l : lOssenkopf 



1993 



Kimura et al. 



2006 



Bertini et al.ll2007l ) 



also known as "Ballistic Particle-Cluster Aggregation" (BPCA). Clusters formed in this way have 
fractal dimension ~ 3, but high porosity. Figure [T] (top row) shows examples of BA clusters with 
N = 256, 1024, 4096 monomers. 

BAMl ("ballistic agglomeration with one migration") clusters are produced by requiring ar- 
riving monomers j > 3, after making first contact with a monomer k < j, to "migrate" to make 
contact with another monomer, by rolling or sliding over the first-contacted monomer, along the 
shortest possible trajectory. If there is more than one candidate for this second contact, the nearer 
is chosen. BAMl clusters with N > 3 have every monomer in contact with at least two other 
monomers; for > 4 some of the monomers are in contact with three or more other monomers. 
Figure [1] (middle row) shows BAMl clusters with = 256, 1024, 4096 monomers. 

BAM2 clusters are constructed as follows: monomers j = 2 and 3 are added randomly just 
as for the BAMl clusters. Monomers j > 4 arrive on random rectilinear trajectories; after first 
contact they make two migrations. The first migration is the same as for constructing BAMl 
clusters: "rolling" along the shortest possible trajectory to make a second contact. This is followed 
by a second migration, now rolling over both the first and second sphere contacted to contact a 
third neighbor, again choosing the shortest trajectory if there is more than one candidate. BAM2 
clusters with > 4 have every monomer in contact with at least three neighbors, with some in 
contact with four or more neighbors. Figure [1] (bottom row) shows examples of BAM2 clusters 
with N = 256, 1024, 4096 monomers. 
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For clusters consisting of N spherical monomers of radius we have, 

aeff = iV^/^ao . (17) 

We generated clusters using the BA, BAMl, and BAM2 algorithms, for a large number of 
different seeds for the random number generator. A library of samples is availabl^ for A'^ = 
2^, 2^, 2^^. Below we report the statistical properties of such clusters. 



3.2. Size and Porosity of Ballistic Aggregates 

The clusters formed following the above procedure are generally irregular and "porous" , with 
the porosity decreasing from BA to BAMl to BAM2, as is apparent from comparison of the three 
rows in Fig. [H 

For a given N and aggregation rule, the characteristic size Rate will vary from realization 
to realization because the aggregation process is random. Figure [2] shows the expectation value 
(Rabc) based on many random realizations, for N = 2^, 2^, 2^^; the realization-to-realization 
variations are also shown. The BA clusters have (Rate) ~^ 2.0aefT as — > oo, the BAMl clusters 
have (Rabc) l-7aeff, and the BAM2 clusters have {Rate) l-5acfr- Note, however, that even for 
N > 10^, the {Rabc/o-cs) continues to increase as N increases from 2^^ to 2^^ to 2^^, so the limiting 
values for N ^ oo are uncertain. Figure [2] also shows {Rkbm), the characteristic size given by eq. 
(fT3l) . For > 10^, i?KBM ~ Rabc- However, for smaller A^ -Rkbm is noticeably larger than Rate- In 
addition, the realization-to-realization variation in -Rkbm is appreciably larger than for Rabc- We 
consider Rabc to be the better way to characterize the effective size of a random structure. 

For each of our clusters, we have estimated the orientation-averaged projected area, and from 
this the radius -Rproj- These computations are time-consuming; the computation is carried out by 
enclosing the target within a sphere of radius Rc, choosing Nr points at random on the sphere, for 
each such point choosing a random inward direction (drawn from a distribution corresponding to 
isotropic incidence on the sphere) , and determining whether that ray does or does not intersect the 
target structure. If the number of intersecting rays is -A^^;, then -Rproj = \/ N^/N^Rc- We typically 
employ = 10^ in order to achieve accuracy better than 1% in the determination of -Rproj- It is 
striking that, for a given .A^, -Rproj in Figure [2] shows less realization-to-realization veriation than 
even Rabc- 

The expectation values (/) and (P) = 1 — (/) are shown in Figure El For large -A^, the BA 
clusters have (P) ~ 0.87, whereas the BAM2 clusters have substantially lower porosity, (P) ~ 0.70. 

Validation of a code to generate random clusters is not simple, but can be done by comparing 
the statistical properties of the clusters with those of other Monte Carlo cluster generators. Statis- 



Ihttp: //www. astro .princeton. edu/$\sim$draine/agglom. html 



Table 1. Sizes of BA, BAMl, and BAM2 Clusters 



Characteristic Radius Rabc/dcH " from eq. (fTTjl Radius -Rkbm " from eq. (fTH)l Projected Area Radius 7i, 



iV BA BAMl BAM2 BA BAMl BAM2 BA BAMl BAM2 



8 


1, 


.4300 


± 


0. 


.0732 


1, 


.2435 


± 


0. 


.0287 


1, 


.1569 


± 


0, 


.0136 


1 


.7049 


± 


0, 


.1683 


1 


.3786 


± 


0, 


.0650 


1 


.2133 


± 


0, 


.0230 


1 


.2569 


± 


0, 


.0163 


1 


.1979 


± 


0, 


.0116 


1 


.1458 


± 


0, 


.0057 


16 


1, 


.5885 


± 


0. 


.0868 


1, 


.3304 


± 


0. 


.0367 


1, 


.1860 


± 


0, 


.0167 


1 


.8209 


± 


0, 


.1671 


1 


.4611 


± 


0, 


.0828 


1 


.2526 


± 


0, 


.0453 


1 


.3514 


± 


0, 


.0197 


1 


.2670 


± 


0, 


.0158 


1 


.1861 


± 


0, 


.0106 


32 


1, 


.7060 


± 


0. 


.0884 


1, 


.4079 


± 


0. 


.0440 


1, 


.2212 


± 


0, 


.0228 


1 


.8906 


± 


0, 


.1486 


1 


.5242 


± 


0, 


.0838 


1 


.2910 


± 


0, 


.0521 


1 


.4469 


± 


0, 


.0219 


1 


.3362 


± 


0, 


.0183 


1 


.2288 


± 


0, 


.0138 


64 


1, 


.7861 


± 


0. 


.0794 


1, 


.4733 


± 


0. 


.0454 


1, 


.2614 


± 


0, 


.0271 


1 


.9262 


± 


0, 


.1227 


1 


.5684 


± 


0, 


.0751 


1 


.3258 


± 


0, 


.0501 


1 


.5405 


± 


0, 


.0227 


1 


.4036 


± 


0, 


.0194 


1 


.2727 


± 


0, 


.0155 


128 


1, 


.8371 


± 


0. 


.0668 


1, 


.5262 


± 


0. 


.0427 


1, 


.2998 


± 


0, 


.0284 


1 


.9403 


± 


0, 


.0966 


1 


.5991 


± 


0, 


.0628 


1 


.3530 


± 


0, 


.0450 


1 


.6299 


± 


0, 


.0225 


1 


.4681 


± 


0, 


.0193 


1 


.3146 


± 


0, 


.0165 


256 


1, 


.8705 


± 


0. 


.0535 


1, 


.5650 


± 


0. 


.0356 


1, 


.3340 


± 


0, 


.0271 


1 


.9464 


± 


0, 


.0745 


1 


.6190 


± 


0, 


.0498 


1 


.3759 


± 


0, 


.0381 


1 


.7138 


± 


0, 


.0216 


1 


.5277 


± 


0, 


.0182 


1 


.3552 


± 


0, 


.0165 


512 


1, 


.8914 


± 


0. 


.0418 


1, 


.5932 


± 


0. 


.0296 


1, 


.3633 


± 


0, 


.0245 


1 


.9470 


± 


0, 


.0565 


1 


.6325 


± 


0, 


.0396 


1 


.3949 


± 


0, 


.0323 


1 


.7894 


± 


0, 


.0199 


1 


.5814 


± 


0, 


.0171 


1 


.3940 


± 


0, 


.0161 


1024 


1, 


.9052 


± 


0. 


.0329 


1, 


.6157 


± 


0. 


.0243 


1, 


.3899 


± 


0, 


.0202 


1 


.9454 


± 


0, 


.0429 


1 


.6428 


± 


0, 


.0297 


1 


.4130 


± 


0, 


.0254 


1 


.8556 


± 


0, 


.0183 


1 


.6289 


± 


0, 


.0158 


1 


.4315 


± 


0, 


.0144 


2048 


1, 


.9177 


± 


0. 


.0248 


1, 


.6343 


± 


0. 


.0183 


1, 


.4114 


± 


0, 


.0171 


1 


.9467 


± 


0, 


.0322 


1 


.6535 


± 


0, 


.0225 


1 


.4274 


± 


0, 


.0207 


1 


.9131 


± 


0, 


.0160 


1 


.6708 


± 


0, 


.0141 


1 


.4647 


± 


0, 


.0130 


4096 


1, 


.9286 


± 


0. 


.0198 


1, 


.6500 


± 


0. 


.0151 


1, 


.4300 


± 


0, 


.0135 


1 


.9493 


± 


0, 


.0250 


1 


.6630 


± 


0, 


.0175 


1 


.4409 


± 


0, 


.0159 


1 


.9612 


± 


0, 


.0148 


1 


.7060 


± 


0, 


.0126 


1 


.4934 


± 


0, 


.0110 


8192 


1, 


.9384 


± 


0. 


.0151 


1, 


.6631 


± 


0. 


.0115 


1, 


.4462 


± 


0, 


.0114 


1 


.9532 


± 


0, 


.0174 


1 


.6718 


± 


0, 


.0130 


1 


.4533 


± 


0, 


.0122 


2 


.0002 


± 


0, 


.0130 


1 


.7344 


± 


0, 


.0101 


1 


.5177 


± 


0, 


.0095 


16384 


1, 


.9508 


± 


0. 


.0121 


1, 


.6758 


± 


0. 


.0097 


1, 


.4588 


± 


0, 


.0082 


1, 


.9606 


± 


0, 


.0129 


1, 


.6828 


± 


0, 


.0113 


1, 


.4640 


± 


0, 


.0088 


2 


.0334 


± 


0, 


.0110 


1 


.7574 


± 


0, 


.0093 


1 1 


.5363 


± 


0, 


.0068 


32768 


1, 


.9620 


± 


0. 


.0089 


1, 


.6875 


± 


0. 


.0073 


1, 


.4714 


± 


0, 


.0063 


1 


.9691 


± 


0, 


.0091 


1 


.6926 


± 


0, 


.0089 


1 


.4751 


± 


0, 


.0066 


2 


.0586 


± 


0, 


.0079 


1 


.7762 


± 


0, 


.0073 




.5531 


± 


0, 


.0065 


65536 


1, 


.9771 


± 


0. 


.0076 


1, 


.6996 


± 


0. 


.0056 


1, 


.4835 


± 


0, 


.0036 


1 


.9811 


± 


0, 


.0089 


1 


.7032 


± 


0, 


.0067 


1 


.4862 


± 


0, 


.0040 


2 


.0795 


± 


0, 


.0086 


1 


.7924 


± 


0, 


.0063 


'l 


.5664 


± 


0, 


.0041 



" Mean and zblcr realization-to-realization variation in Rabc from eq. (|ll|l . 
" Mean and ilcr realization-to-realization variation in -Rkbm from eq. (|13p 
* Mean and ilcr realization-to-realization variation in -Rproj from eq. (|15|l . 
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tical properties of the three types of aggregate clusters are given in Tables [T] and [2] for cluster sizes 
ranging from = 2^ to = 2^^. 



The BA clusters have been frequently used in th e literature (e.g., iKozasa et al 



1992 



1993 



Ossenkopi 


1993; 


Kimura et al. 


2006: 


Bertini et al. 


2007) 



areas agree with results for BA clusters reported by iKozasa et al.l (|l993l ) and iNakamura et al 



(|l994l ) as shown in Figure [2] (top panel). However, as seen in Figure [2](mid dle panel), ou r value s 
of -Rkbm (defined by eq. [T3|) are slightly lower than the results reported by iKozasa et al.l (|l993l ) . 
The reason for this discrepancy is not known I§ 



3.3. Shape of Ballistic Aggregates 

The axial ratios c/h and h/a give an indication of the overall shape of each random aggregate. 
Figure SI shows c/h and h/a for BA, BAMl, and BAM2 clusters with A^ = 2^, 2*^, 2^0, and 2^2 
spheres. For a given A^, the sequence BA— >BAM1— >BAM2 corresponds to decreasing eccentricity. 
Similarly, for a given agglomeration rule (BA, BAMl, or BAM2), clusters with larger A^ tend to 
be more spherical. 



T. Kozasa has kindly provided several BA clusters from IKozasa et al. For these clusters, we obtain 



values of Rabc/cLcH and -RKBM/ffloff that are consistent with the statistics reported in our Table [T] 
Table 2. Porosities of BA, BAMl, and BAM2 Clusters 



Porosity V " from eq. (|12|l Porosity Pkbm 



TV BA BAMl BAM2 BA BAMl BAM2 



8 





6580 


± 





0525 





4799 


± 





0360 





3541 


± 





0228 





7982 


± 





0598 





6183 


± 





0540 





4402 


± 





0319 


16 





7505 


± 





0409 





5753 


± 





0351 





4006 


± 





0253 





8265 


± 





0443 





6736 


± 





0519 





4874 


± 





0522 


32 





7986 


± 





0313 





6417 


± 





0336 





4509 


± 





0307 





8468 


± 





0338 





7127 


± 





0444 





5309 


± 





0532 


64 





8245 


± 





0234 





6873 


± 





0289 





5018 


± 





0321 





8568 


± 





0259 





7374 


± 





0359 





5674 


± 





0467 


128 





8387 


± 





0176 





7187 


± 





0236 





5446 


± 





0299 





8611 


± 





0200 





7533 


± 





0280 





5937 


± 





0390 


256 





8472 


± 





0131 





7391 


± 





0178 





5788 


± 





0257 





8632 


± 





0152 





7630 


± 





0213 





6144 


± 





0312 


512 





8522 


± 





0098 





7527 


± 





0138 





6053 


± 





0213 





8638 


± 





0115 





7694 


± 





0164 





6304 


± 





0250 


1024 





8554 


± 





0075 





7629 


± 





0107 





6276 


± 





0162 





8638 


± 





0089 





7740 


± 





0121 





6449 


± 





0189 


2048 





8582 


± 





0055 





7709 


± 





0077 





6443 


± 





0129 





8642 


± 





0066 





7786 


± 





0090 





6558 


± 





0149 


4096 





8606 


± 





0043 





7774 


± 





0061 





6580 


± 





0097 





8649 


± 





0051 





7824 


± 





0068 





6655 


± 





0110 


8192 





8627 


± 





0032 





7826 


± 





0045 





6694 


± 





0078 





8657 


± 





0036 





7859 


± 





0050 





6470 


± 





0082 


16384 





8653 


± 





0025 





7875 


± 





0037 





6779 


± 





0054 





8673 


± 





0026 





7901 


± 





0042 





6812 


± 





0057 


32768 





8676 


± 





0018 





7919 


± 





0027 





6861 


± 





0040 





8690 


± 





0018 





7937 


± 





0032 





6884 


± 





0042 


65536 





8706 


± 





0015 





7963 


± 





0020 





6937 


± 





0022 





8714 


± 





0017 





7976 


± 





0024 





6953 


± 





0025 



" Mean and ±1(t realization-to-realization variation in porosity V from eq. (|12|l . 
^ Mean and ±lcr realization-to-realization variation in porosity Pkbm from eq. (|14l ') 
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Prolate spheroids would have h/a = 1; oblate spheroids have c/h = 1. Figure U] shows that the 
clusters are generally triaxial, with a tendency toward prolateness {c/h > b/a). 

The ballistic aggregates constructed in this section are based on very simple algo r ithms. More 



realis tic treatments of grain coagulation processes are possible (e.g., IOssenkopflll993l : lOrmel et al 



20071 ). but are rather uncertain given our limited knowledge of the interstellar /circumstellar envi- 
ronment and grain properties. The well-defined sequence of BA to BAM2 clusters described here 
allows us to systematically investigate the grain properties as function of porosity. 



4. Application of the Discrete Dipole Approximation 



4.1. Composition 

The composition of interstellar grains continues to be uncertain. The observed strength of the 
10 fim absorption feature, produced by the Si-0 stretching mode in silicates, requires that >50% of 
the total grain volume be amorphous silicate material. As discussed below, this amount of silicate 
material is also consistent with the observed depletions of Si, Mg, and Fe from the gas phase. 

Attempts to reproduce the observed wavelen gth-dependent extinction require more grain ma - 



terial than can be provided by silicates alone (e.g, IWeingartner &: Draing 



20011: 



Zubko et al.ll2004l '). 



Because H2O ice is not present in the diffuse interstellar medium, and Mg, Fe, and Si are presumed 
to be primarily invested in the silicate material, the only element that can provide substantial 
additional sol id material is carbon. Observational evidence for carbonaceous material includes (see 
Draind l2003al . and references therein): (1) strong absorption near 0.22 //m, likely due to vr — > vr* 
electronic excitation in sp^-bonded carbon (e.g., graphite, or polycyclic aromatic hydrocarbons); 
(2) absorption at 3.4 /im (attributed to the C-H stretching mode); and (3) observed emission fea- 
tures at 3.3, 6.2, 7.6, 8.6, 11.3, and 12.7 ^m that are attributed to vibrational modes of polycyclic 
aromatic hydrocarbons (PAHs). The cosmic abundance of carbon, and the fact that it is mod- 
erately depleted from the gas phase, is consistent with carbonaceous material with a total solid 
volume that is a substantial fraction of the volume of amorphous silicate material. 

Assuming total interstellar abundances of the elements relative to H to be the same as current 
estimates of solar abundances, observations of gas-phase abundances in the inters tellar inedium 
allow us to infer the amount of different elements locked up in grains (see Table 1 of iDraind l2008l ) . 
Current estimates indicate that carbon in grains contributes a mass of about 0.13% of the total 
H mass, while solid material containing Mg, Fe, Si, O (plus small contributions from Al, Ca, and 
Ni) amounts to about 0.57% of the total H mass. If the carbonaceous material has a density of 
~ 2gcm~'^, and the silicate material has a density ~ 3.5 gcm^'^ then the overall silicate/carbon 
volume ratio is Ku/^car ~ 2.5. Thi s, however, is based o n assuming the total C abundance to 
be only 245 ifc 30 ppm relative to H (lAsplund et al.ll2005bl). H owever, two recent determinations 
of (O/H)0 (ILandi et al.l 120071 ; ICenteno &: Socas-Navarrd 120081 ) obtain values that are ~ 1.9 times 
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larger than the solar oxygen abundance (O/H)© = 457 ± 56 ppm of lAsplund et al.l (j2004l . l2005bl ). 
The solar C abundance might therefore be larger than Asplund et al.'s value of 245 ppm. If the total 
C abundance were to be increased to, e.g., 350 ppm, the mass of C in dust would increase by 75%, 
and the silicate/carbonaceous volume ratio would fall to Vsn/Vcur ~ 1-4. The total abundances of 
Mg, Fe, Si - and therefore the inferred abundance of solids with silicate composition - are of course 
also uncertain. We conclude that if grains are of mixed composition, the silicatexarbonaceous 
volume fractions could be as large as 70:30, or as low as 60:40. We will consider extreme examples 
of composite grains that are 100% silicate, as well as composite grains with 50:50 volume fractions. 

While the importance of carbon is undisputed, the specific form is uncertain. The observed 
PAH emission features require ^20% of the solid carbon to be in small free-flying PAH molecules 
or clusters. Spectroscopy of the 3.4 fim feature indicates th at ^ 15% of the hydrocarbon ni aterial 
is aliphatic (chain- like) , and ~ 85% aromatic (sp'^-bonded) ( Pendleton Allamandolall2002l ). Dia- 
mond (sp^-bonded carbon) has been found in meteorites, and may also be present in the interstellar 
medium, but there is no direct evidence for interstellar diamond. 

Because the predominant form of carbon is uncertain, we consider two possibilities: crystalline 
graphite, ar id amorphous c arbon. For graphite, a highly anisotropic material, we use the dielectric 
tensor from lDraind (j2003bl ). Each graphite sphere is assumed to be a single crystal, with the crystal 
axes for each graphite monomer assigned an independent random orientation. 

"Amorphous carbon" is not a well-defined material, and its properties depe nd on the method of 
preparation. We use the dielectric function for amorphous carbon "ACl" from lRouleau &: Martin 
(1199 Ih. For th e amorphous silicate component, we use the "astrosilicate" dielectric function from 
Drainel (|2003bl ). 



We will consider three compositions for the clusters: (1) 50% amorphous silicate and 50% 
graphite; (2) 50% amorphous silicate and 50% ACl; and (3) 100% amorphous silicate. For the 
mixed-composition clusters, the monomer compositions are assigned randomly, and the percentages 
are by volume. 



4.2. DDA Method and Validity Criteria 



We calculate th e absorption and scatte ring properties of the BA, BAMl and BAM2 clusters 



proximation (jPurcell Pennypacker 



197: 



I 



using DDSC AT 7.0 (IDraine fc Flataull200^ ). DDSCAT is a code based on the discrete dipole ap- 



Draine &: Flatau]ll994l ). designed to compute scattering 



and absorption of electromagnetic waves by targets with arbitrary geometry and composition, for 
targets that are not too large compared to the wavelength A. 



There are three validity criteria that should be satisfied for the DDA to provide accurate 
results: 
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1. The interdipole spacing d should be small enough to resolve the geometric structure of the 
target. This is accomplished provided the number of dipoles per spherical monomer n^ip ^ 
100. 

2. The phase shift corresponding to one dipole spacing should be small. 

MM = 044|»| f ™ ) (^^) ("-Lt^) < 1 (18) 

where m is the refractive index, k = 27r/A, ao is the monomer sphere radius, and n^ip is 
the number of dipoles per sphere. Thus we see that n^ip = 100 allows condition ([18]) to be 
satisfied even at A = 0.1 jim^ for oq = 0.02 /im. For accurate calculations of the scattering 
phase function, the validity requirement is somewhat more stringent, \m\kd < 0.5. With 
^dip = 100 and ao = 0.02 /im, we have \m\kd < 0.5 for A > 0.2 /xm. 

3. Even when the criterion \m\kd < 0.5 is satisfied, the implementation of the DDA used here 
overestimates absorption in materials with Im(e) ^ 1, where e is the dielectric function. 
Graphite is a conducting material; the dielectric tensor component for E J- c (with c being the 
normal to the basal plane) becomes large in the infrared (e_L = 11.7 + 32. 6z at A = 3.981 /um), 
and the accuracy of the DDA suffers. 



To assess the computational accuracy, some of the scattering calculations have been repeated using 
differing numbers of dipoles per monomer n^ip. In the limit n^ip oo, the DDA becomes exact. 

— 1/3 

For finite njip, the errors are expected to scale as n^^^ (i.e., the errors var y linearly with the 



interd ipole separation d) ; such behavior has previously been demonstrated bv ICollinge &: Draine 



\200i ] and lYurkin et al.l ((20061). 



Figure [5] shows the total extinction and scattering efficiencies computed for one cluster (50% 
silicate and 50% graphite) with N = 256 ao = 0.02 /xm monomers, and porosity V = 0.853, for a 
single orientation, using different numbers of dipoles with njip ~ 100, 200, 400, 1000. The results 
show that the computed cross sections at each wavelength are approximately linear functions of 
n^^J^, allowing us to confidently extrapolate to estimate the exact result at n^^J^ — s- 0, and 
thereby to estimate the error that would result if we were to simply use cross sections calculated 
for ndip ~ 100. On average, the fractional errors for the results with n^ip ~ 100 are only a few 
percent for A < 0.35//m, 4-10% for 0.35 < A < 0.8/im, and 10-15% for 0.8 < A < 4//m. Note that 
for wavelengths A > 0.5 /im, the computations for finite ndip always overestimate the absorption 
cross section and scattering cross sections, presumably as a result of failure to accurately resolve 
the "shielding" produced by the charge layer produced by the discontinuity in the polarization field 
at the grain surface. 



''At A — 0.1 /im, |m| = 1.84, 2.41, 2.15, 1.91 for astrosilicate, graphite with E ± c, graphite with E \\ c, and ACl. 
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The fractional errors of ~ 10% in the infrared are comparable to the typical variations in Qabs 
and Qext from one random realization to another. Hence using njip ~ 100 dipoles per monomer 
suffices for most of our investigation^, and allows us to explore the parameter space more efficiently. 

For the DDA results presented in the next section, each calculation is averaged over a few (3 
or 5) realizations and 54 random orientations for each realization. The typical variations in cross 
sections from realization to realization for the same type of clusters are of the order < 10%; while 



54 ra ndom orientations are sufficient to represent the orientational average (see I Johnson &: Draine 



20081 ). The scattering and absorption are calculated at 32 wavelengths from 0.1 to 3.981 //m (with 
AlogioA = 0.05, 0.025, 0.05, 0.1 for 0.1 < A < 0.1259, 0.1259 < A < 0.3126, 0.3126 < A < 0.5012 
and 0.5012 < A < 3.981 fim, respectively). 



5. Absorption, Scattering, and Extinction Cross Sections of Aggregates 

In § §5.1|5.2l we fix the effective radius a^s = 0.127 fim for each cluster and explore the differ- 
ences in optical properties between three different cluster geometries and three different composi- 
tions. The aggregates have characteristic radius Rate ~ 0.16— 0.24;^m. In §5.2l we vary the monomer 
size to investigate the sensitivity of results to monomer size and to porosit y at fixed Cpff. In i^5.3lwe 



comp are our DDA results with the analytical multi-layer sphere model (jVoshchinnikov &: Mathii 



np 

3) 



19991 ) to test the accuracy of the MLS prescription for estimating the optical properties of ballistic 



aggregates. In ^5.41 we carry out a similar comparison of our DDA results with "effective medium 
theory" . 



5.1. DDA Results: Dependence on Wavelength and Porosity 

The orientation- and realization-averaged (5 realizations) extinction cross section, absorption 
cross section and scattering cross section, as well as the asymmetry parameter g = (cos 6) are plotted 
in Figure [6l for the three aggregate types (BA, BAMl, BAM2) and three different compositions, 
as functions of wavelength. We have used the results computed with n^ip ~ 400, which are within 
a few percent of the exact values as inferred from Fig. [5l 

For wavelength A ^ a^s, the clusters are in the Rayleigh limit; the extinction is dominated 
by absorption, and the asymmetry parameter g = (cos 9) is small. For the same volume of solid 
material, the 100% silicate clusters have smaller absorption (and therefore total extinction) at long 
wavelength, because graphite and ACl are more absorptive than silicate in the near-infrared. At 
short wavelength, the behavior of the cross sections and g is nonmonotonic. For example, for the 
50% silicate/50% graphite clusters, Qext bas a dip at ~ 0.17/im and peak near 0.22 /xm. 



^ For clusters with A*' < 512, we use results obtained with ridip ~ 400 for improved accuracy; for larger clusters, 
i.e., A'^ — 1024, 2048, etc., we use ridip ~ 100 due to computational limits. 
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An important question is how the porosity of the dust grain affects the extinction efficiency. 
One might expect that increased porosity would lead to an increase in the overall extinction cross 
section per unit solid material, which motivated attempts to try to use porous dust grains to 
circumvent the difficulties in accounting for the observed interstellar extinction without overcon- 
suming the elements used to build grains. Fig. [6] shows that for short wavelengths, increased 
porosity does result in a modest increase in extinction. However, this effect reverses longward of a 
transition wavelength Xf. The transition wavelength Xt ~ 0.34 ^um for 50% silicate/50% graphite, 
~ 0.25 ^m for 50% silicate/50% ACl, and ~ 0.20 i-im for 100% silicate. For these three examples, 
Xt « (2.1±0.6)aeff. 

Although in this section we focused on clusters with = 0.127//m, similar behavior is seen 
for other cluster sizes. The transition wavelength at which the effect of porosity reverses increases 
when the overall size of aggregate clusters increases. Based on limited numerical experiments, we 
conjecture that, in general, the transition wavelength occurs at Xt ~ (2.5±1.0)oefr: porosity appears 
to increase the extinction cross section per unit solid material for A < l.Saeg, while reducing the 
extinction per unit solid material at A > 3.5acfT. 

It is interesting to no te that for the 50% silicate/50% ACl case, similar trends are seen in 



Voshchinnikov et al. 



(120061 . fig. 2), where the porous dust grain is spherical and modeled using the 
EMT-Mie theory. Th e fact t hat p orosity actually decreases the extinction efficiency at A > A^ ~ 
2aefr was also seen by IWestI ()l99ll ) , who used two types of aggregates with different porosity and 
compared with the results of equal-volume spheres, and found that porosity caused a decrease in 
extinction efficiency when the equivalent size parameter is IttR/X < 5 (i.e., fig. 3 in Iwestl[l99ll ). 



5.2. Effects of monomer size/sensitivity to porosity 

In the above discussion, we fixed the monomer size to be 200 A in radius. There is no direct 
evidence of how large the monomers should be, hence we need to know if different monomer size 
inside a ballistic aggregate will affect the overall absorption and scattering cross sections. To this 
end, we choose two different monomer sizes, 100 A and 400 A, for the BAM2 clusters with 2048 
monomers and 32 monomers respectively. These results will be compared with our fiducial BAM2 
clusters with 256 ao = 200 A monomers discussed in ^5. It we are comparing three clusters with the 
same Oeff, i.e., composed of the same amount of solid materials. We focus on the 50% silicate/50% 
graphite case. For every cluster we average over three realizations, and many random orientations 
of each realization. We used n^ip ~ 100 for the 2048-monomer clusters and n^ip ~ 400 for the 
256-monomer and 32-monomer clusters. To satisfy the DDA accuracy criterion equation (jlSp we 
only focus on results with wavelength A > 0.2 /im. We note that although the three clusters contain 
identical volumes of solid materials with aeg = 0.127 //m, the porosity is different for the 32, 256, 
2048 clusters, with V k, 0.45, 0.58, and 0.64, respectively. 

We show the results in Fig. [3 for Qcxti Qabs) Qsca and g = {cos 6) respectively. Since the 
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three clusters have same amount of sohd materials but different porosity (porosity increases with 
increasing number of monomers), this behavior is very similar to the general behavior of BA, BAMl 
and BAM2 clusters as shown in Fig. [6l Once again, there is a "transition" wavelength At ~ Soefr 
where the calculated Qext is approximately independent of V (the Oefr = 0.127 ^m clusters in Fig. [6] 
have At ~ 0.34 //m 2.7aeff, and the a^E = 0.127 /xm clusters in Fig.[7]have At ~ 0.42 /^m S.Saefj). 

Although in Fig. [7] there seems to be a second transition wavelength at ~ 1.5 /im, we cau- 
tion that it might be artificial because the accuracy of our DDA results decreases at the longest 
wavelengths (see Fig. [5|) for small n^ip. This is particularly true for the 2048 case where we used 
™dip ~ 100 due to computational limits. 

To further demonstrate that the effect seen in Fig. [7] is actually the result of varying the 
porosity, we carry out another test. In Figure [8] we compare two clusters with different geometries 
(BAMl vs. BAM2), and different monomer sizes (504 A vs. 200 A), but with the same amount 
of material (aeg = 0.160;um), and approximately the same porosity {V ~ 0.62). The two clusters 
have very similar cross sections for extinction, absorption, and scattering. Thus we conclude that 
the important parameters are just a^s (i-e-, the amount of solid material) and the porosity 7-", and 
monomer size does not have significant effects as long as the monomer size oq ^ A/27r. 



5.3. Comparison with the MLS Approximation 

The DDA computations described above are very time-consuming. To model p orous, composite 



dust grains, some authors have used th e analytical multi-layer sphere (MLS) model (jVoshchinnikov Mathis 



199S 



Voshchinnikov et al.ll2005l . 120061 ). where the sphere is composed of concentric spherical shells, 
each of which is further composed of a set of spherical layers of single composition. The number of 
such shells should be large enough so that the results are unaffected by changing the order of layers 
inside each shell; in other words, the materials are well-mixed inside such spheres. The problem 
of light scatteri ng and absorption by a multi-layer sphere can be solved by a fast algorithm (e.g.. 



Wu et al.l 119971 ). The applicability of the MLS model to the aggregates considered here has not 



previously been examined. 

To compare the DDA results with the MLS model, one needs to know the vacuum fraction /vao 
which is the volume fraction of vacuum in the multi-layer sphere. Proponents of the MLS method 
have not addressed the question of what value to use for /vac when modeling random aggregates; 
the optimal value of /vac is not necessarily equal to the porosity V given by eq. (fT2]) . Therefore we 
test the MLS using different values of /vac to find the value of /vac that minimizes the difference 
between the DDA results and the MLS results. In practice we set /vac = 0.1 — 0.9 with an increment 
of A/vac = 0.1 in the MLS calculations. We set the number of shells in the MLS to be 40 so that 
the results do not vary significantly when the order of layers inside each shell is changed. 

We plot the results of MLS calculations on top of the DDA results in Fig. [6] as black lines, which 
are bounded by the black dashed line with /vac = 0.1 and the black dotted line with /vac = 0.9. MLS 
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results using intermediate vacuum fractions lie between these two bounds. One can immediately 
see that the MLS model is not a good approximation, especially for the 50% silicate/50% graphite 
case. For the 50% silicate/50% ACl case, the MLS results have two transition wavelengths (~ 
0.31 fim and ~ 0.95 am) a t which the effect of porosity reverses, as has already been reported by 



Voshchinnikov et al.l (|2006l ). although Oeff = 0.1 fim in their work and 0.127/im here. However, this 
is different from our DDA results where there is only one transition wavelength, near ~ 0.25 /xm. 
Similar results are found for the 100% silicate case. For the 50% silicate/50% graphite case the 
MLS model fails to predict such a transition altogether - increasing the porosity increases the MLS 
extinction at all wavelengths between 0.1 and 4//m. For all three compositions, the MLS model 
seems to underestimate the extinction at short wavelengths and overestimate the extinction at long 
wavelengths. Our DDA results are of limited accuracy (errors of up to ~ 8%) at long wavelengths 
for n^ip ~ 400 (see Fig. [5]), but the exact extinction values are even lower at these wavelengths, 
hence the MLS errors are even larger than shown in Fig. [6l 

To quantify the deviation of the MLS results from the DDA results we define a global error 
averaged over all 32 wavelengths: 

(global error Q^f = {[ln{Q^^MLs/Qx,DDA)?) ■ (19) 

Fig. [9] shows the global errors of Qcxt, Qabs, and Qsca for the three types of aggregates and 
three compositions, averaged over 5 realizations. Dotted lines show the standard deviation from 
the 5 realizations. It is clear that the accuracy of the MLS model decreases when the porosity 
increases (from BAM2 to BA clusters) or when the constituent materials become more absorptive 
(from 100% silicate to 50% silicate/50% graphite). 



5.4. Comparison with the EMT Approximation 

One approach to estimate the optical properties of random aggregates is to approximate them 
by homogeneous spheres, with an "effective" refractive index obtained from "effective medium the- 
ory" (EMT) ; the scattering and absorption by the homogeneous sphere is then calculated using Mie 
theory. Effective medium theo ry comes in more than one variant; here we consider the form of EMT 



developed by Bruggeman (see lBohren &: Huffmanlll983l ) where the effective dielectric permittivity 
eeff is calculated via 

y /.ii^i£^ = , (20) 

'I 

where /j and ej are the volume fraction and dielectric permittivity of each composition, including 
vacuum. Note that eeff calculated in this way is not affected by the detailed structure of the 
composite grain, i.e., "monomer size" does not matter. If the aggregate contains n distinct materials 
(including vacuum), then eeff is a root of an n*'^ order complex polynomial equation. We find that 
there is always just one root with Im(ecff) > - this is the physically meaningful solution. 
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The cross sections calculated using EMT are compared to the DDA results in Fig. [10] using our 
fiducial N = 256 clusters with uq = 0.02 //m, and the global errors from the EMT approximation 
are shown in Fig. [TTl We take the vacuum fraction in the EMT-Mie model to be /vac = 0.1 — 0.9 
with A/vac = 0.05. The EMT-Mie approach provides much better agreements with the DDA 
results than the MLS approximation does. This is particularly true for the silicate-graphite and 
the silicate-ACl compositions, where one material is highly absorptive. It is clear from Figs. [6] 
and [10] that the dif ference between MLS and E MT-Mie increases when /vac increases, which is also 



evident in fig. 2 of 1 Voshchinnikov et al.l (|2006l ) for the silicate-ACl case. Our ballistic aggregates 



have porosities V > 0.6, hence it is not surprising that the MLS model is not a good approximation. 

The global error plots in Fig. [11] show that the EMT-Mie model gives the optimal results 
wheM /vac = 0.80, 0.70, 0.55 for BA, BAMl, BAM2 clusters respectively: it appears that for the 
three types of aggregates the optimal vacuum fraction for EMT-Mie calculations is 

/vac ~ 0.94P . (21) 

The global errors of Qcxt etc. are < 10% at this optimal /vac, which suggests that the EMT- 
Mie model provides a fairly good approximation for computing overall cross sections. We plot 
the difference between the EMT-Mie results (for optimal choice of /vac) and the DDA results, as 
functions of wavelength, in Fig. [12] It shows that although the EMT-Mie results do not follow the 
DDA results exactly, the maximum deviations are typically ~ 20% in all cases. It is also apparent 
that the EMT-Mie model produces more forward scattering at short wavelength, a property that 
is further discussed in Paper II, which examines the angular distribution and polarization of the 
scattered light. 



6. Summary and Discussions 

The principal results of this study are the following: 

1. Two new algorithms for generating random aggregates are introduced: ballistic aggregation 
with one migration (BAMl), and ballistic aggregation with two migrations (BAM2). BAMl 
and BAM2 aggregates are less porous, and more mechanically robust, than conventional BA 
aggregates. 

2. A measure V of the porosity of a structure is proposed (eq. [T0|) . as well as a measure Rabc 
(eq. [TT|) for the "characteristic size" of the structure. 

3. Monte-Carlo simulations are used to determine the statistical properties of V and Rate for 
ballistic aggregates (see Tables 1-2 and Figs. [2][3]). 



These are not the exact values of porosity at which the EMT-Mie model fits the DDA results best, because our 
grid of vacuum fraction has a coarse grid size of A/vac ~ 0.05. 
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4. We confirm (see Figure [5|) that the accuracy of the DDA scales as the interdipole spacing d 
as d — > 0, or, equivalently as N^^^^^ as A'dip oo, where A^dip is the total number of dipoles. 

5. Scattering, absorption, and extinction cross sections are calculated for the N = 256, oq = 
0.02 /um (oeff = 0.127 /im) BA, BAMl, and BAM2 clusters, for three different compositions: 
100% silicate, 50% silicate/50% amorphous carbon, and 50% silicate/50% graphite, for wave- 
lengths 0.1 < A < 3.981 fim. The BA clusters (with the highest porosity V) have the largest 
extinction cross sections at short wavelengths, but at optical and near-IR wavelengths, the 
BAM2 clusters (with the lowest V) provide more extinction per unit solid material. At con- 
stant porosity and same amount of solid material, the monomer size has no significant effect 
provided the monomers are small compared to the incident wavelength. 

6. We compared the DDA results with the analytical MLS model and EMT-Mie theory. We 
found the MLS model does not provide a good approximation for absorptive and/or very 
porous grains; the EMT-Mie model provides much better agreement with the DDA results. 
For computing total cross sections {Qext, Qabs; Qsca), the EMT-Mie method provides results 
accurate to ~ 10% if the vacuum fraction /vac is taken to be 0.947^. 



The effects of porosity on extinction cross sections have important implications for the abun- 
dance budget_2roy^mjnmterstellar dust models. The recent observed decrease of solar abundances 



(e.g.. lAsplund et al.ll2005al) and the claim th at interstellar abundances might be better represented 



by those of B stars (e.g., ISnow Wittlll996l ) have imposed a challenge to dust extinction models. 
Porous dust grains have been thought to be a solution to this abundance budget problem (e.g., 
Mathid 119961 ). because they were expected to result in greater extinction per unit solid material 
than compact grains. However, our results show that porosity actually decreases the opacity at 
wavelengths long compared to the overall grain size. Hence caution must be paid when dealing 
with the abundance budget problem. Until detailed models have been constructed using random 
aggregates to reproduce the observed interstellar extinction (and polarization), we will not know 
if such gra in models will alleviate the interstellar abundance problem. Work on this problem is 
underway (| Johnson &: Draind 120081 ). 



Ballistic aggregates are promising candidates for interstellar and circumstellar dust grains. In 
the companion Paper II, we will discuss the scattering properties of ballistic aggregates and present 
examples that can reproduce the observations of light scattered by dust in debris disks and comets. 



We thank the anonymous referee for helpful comments. This research was supported in part 
by NSF grant AST-0406883. Computations were performed on the Delia and Artemis computer 
clusters at Princeton University. 
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Fig. 1. — From left to right: N=256, 1024, and 4096 clusters. From top to bottom: BA, BAMl, and BAM2 
aggregation rules. Axes ai, 3.2, £13 are the principal axes with the largest, intermediate, and smallest moment of 
inertia. For each cluster we give the characteristic radius Rate, the porosity V, and the dimensionless moment-of- 
inertia eigenvalues at. 
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Fig. 2. — Characteristic radius i?afec, -Rkbm, and the radius -Rproj based on mean projected area, for 
clusters produced by BA, BAMl, and BAM2 aggregation rules with single-size spheres (see text). 
Solid lines connect the mean values. Dotted lines show t he ifclo" dispers i on for random realizations. 
For BA clusters , we sh ow -Rkbm and -Rproj reported by iKozasa et al.l (jl993l ). and -Rproj given by 
Nakamura et al.1 (119941 ). 
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Fig. 3. — Volume filling factor / (left scale) and porosity V (right scale) for clusters produced 
by BA, BAMl, or BAM2 aggregation rules with single-size spheres. Solid lines connect the mean 
values. Dotted lines show the itlo" dispersion for random realizations. BA clusters are the least 
dense, with porosity V ~ 0.87 in the limit N oo. BAM2 clusters are the most dense, with 
V ^ 0.70 in the limit iV ^ oo. 
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Fig. 4.— Axial ratios for BA, BAMI, and BAM2 clusters of N spheres, for N = 64, 256, 1024, and 
4096. BA clusters are more asymmetric than BAMI clusters, which in turn arc more asymmetric 
than BAM2 clusters. Prolate shapes (c/6 > b/a) appear to be slightly favored, and there is an 
overall tendency for clusters to be rounder as N increases. 
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Fig. 5. — Convergence tests for the DDA method. Shown here are the extinction (left) and scat- 
tering (right) efficiencies for one N = 256, ao = 0.02 ^m, 50% sihcate and 50% graphite cluster 
with V = 0.853, computed using n^ip « 100, 200, 400, 1000, for selected wavelengths (as marked 
on top- left of each panel, in units of fan.), and for a single orientation. The filled circles are the 
results, and the dotted lines are extrapolations to n^ip 00 using the n^ip ~ 400 and 1000 results 
(see text). The percentage alongside each data point is the fractional error with respect to the 
extrapolated value at n^ip 00. 



-27- 



50% sil + 50% gra 50% sil + 50% AC1 100% si 




0.1 0.2 0.5 1.0 2.0 0.1 0.2 0.5 1.0 2.0 0.1 0.2 0.5 1.0 2.0 
A (yLim) A (yLim) A (/xm) 

Fig. 6. — Wavelength dependence of orientation-averaged Qext, Qnhs, Qsca. and g = (cos 6*) for three 
compositions and three aggregate types, for the fiducial = 256 and ao = 0.02 //m clusters with 
n^ip ~ 400. Each type of cluster is averaged over 5 realizations. Also plotted are the results of the 
MLS model: the black dashed lines are the MLS results with /vac = 0.1, the black dotted lines are 
for /vac = 0.9, and the black solid lines are for /vac = 0.3, 0.5, 0.7. 
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0.2 0.5 1 2 0.2 0.5 1 2 

A (/xm) A (/xm) 

Fig. 7. — Qext ; Qabs , <3sca and g = (cos 9) , for three B AM2 aggregates with 50% sihcate and 50% 
graphite. These aggregates have the same amount of sohd materials {qcs = 0.127/um) but different 
porosities {V ^ 0.45, 0.58, 0.64) for N = 32, 256, 2048. At a given A, the computed Qabs, Qext, 
and Qsca sue porosity-dependent, except at the "transition radii" At ~ 0.35/im and ^ 1.4/Lim where 
the cross sections are insensitive to porosity. For 0.35 < A < 1.4;um, the cross sections decrease as 
the porosity is increased. 
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Fig. 8. — Comparison of two clusters with different geometries (BAMl vs. BAM2), different 
monomer size but similar porosity V « 0.62 and same Oeff = 0.160 //m. The optical properties are 
very similar. Therefore monomer size has no significant effects as long as monomers are smaller 
than the incident wavelength. 
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Fig. 9. — Global errors (defined in eqn. [T9]) for the MLS results and the DDA results for the fiducial 
aggregates as shown in Fig. [6l Dotted lines show the standard deviation from the 5 realizations. 



- 31 - 



50% sII + 50% gra 50% sil + 50% ACl 100% si 




0.1 0.2 0.5 1.0 2.0 0.1 0.2 0.5 1.0 2.0 0.1 0.2 0.5 1.0 2.0 
A (fxm) X (/i.m) A (/xm) 

Fig. 10. — The same as Fig. [6l only with the MLS results replaced by the EMT-Mie results, for 
/vac from 0.1 to 0.9. 
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Fig. 11. — Global errors (similar definition as in eqn. [T9|) for the EMT-Mie results and the DDA 
results for the fiducial aggregates as shown in Fig. \W[ /vac from eq. ([211) is shown. 
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Fig. 12. — The difference between tlie optimal EMT-Mie results with /vac = 0.947^ and the DDA 
results for our fiducial clusters. Plotted here are (Qemt/Qdda — 1) for cross sections, and (^emt — 
Odda) for g = (cos 6*). The fractional difference is typically < 20% for total cross sections. 



